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POLYNOMIAL STRICTION ALGORITHMS IN THE DETECTION AND DETERMIN- 
ATION OE THE CHARACTERISTICS OF NON - ST AT I ONARY PHENOMENA 

V. P. Yevdokimov and M. Ya. Natenzon 


The task of the detection and measurement of the character- / 2* 

istics of shock waves and explosions in the cosmic plasma, which 
lead to the non-stationariness of behavior of the physical para- 
meters measured on board a space vehicle, requires provision of 
the maximum information content of the experiment in the area of 
the intermittent change ?n the parameters, and assumes a lesser 
information content, or even the absence of measurements, in the 
area of their "calm" behavior. Examined herein is the effectiveness 
of the algorithm of polynomial striction of the data, and the 
selection of an algorithm is carried out, according to the results 
of measurements of the components of the permanent magnetic field 
in communication sequences from the "Venera-9" space vehicle, 
where the shock waves were recorded. It is shown that the most 
effective algorithm, according to the coefficient of striction, 
is the algorithm of the zero order interpolator (INP). Also ex- 
amined is the increase in effectiveness with the introduction of 
a variable threshold of comparison. Calculated herein is the 
root-mean-square deviation of the error of regeneration of the 
data with various thresholds, and qualitative analysis of the 
distortions is carried out, having shown a considerable change 
in the form of the small discontinuities, which may be inter- 
preted as potential interplanetary shock waves. With a selected 
magnitude of the threshold 1* , the magnitude of the coefficie it of 
striction proved to be relatively small, and (without taking into 
account the function information) did not exceed six. 

1. Posing of the Problem /3 

An important task of space physics, being solved intensively 
in recent years, is the problem of detection and measurement of the 

*Numbers in the margin indicate pagination in the foreign text. 


1 


characteristics of the shock waves and explosions in the cosmic 
plasma, which lead to non-stationar inoss in the time of the 
physical parameters measured on board a space vehicle. The object 
of the study, in the total volume of obtained data, is some moments 
of intermittent changes in the measured parameter, and the section 
prior to and after these measurements. The obtaining of a full- 
volume of data with a maximum prescribed frequency of sampling, 
which provides reproduction of the parameter with guaranteed 
accuracy, should be ensured in these sections. The utilization 
of such a high frequency of sampling in sections of "calm" 
behavior of the parameter is inadvisable, since it considerably 
increases the volume of information transmitted from onboard 
the space vehicle. 

There occurs the problem of sampling of the information prior 
to its transmission to earth, which may be solved by two methods: 

1. Utilization of a variable frequency of quantization with 
time, adapted to the behavior of the parameter being studied and 
ensuring a guaranteed accuracy of reproduction an earth. 

2. Detection of non-stationary states of a certain type and 
transmissions to earth of data only in their vicinities. 

The former method is related to the so-called "quasiinvertible" 
methods of striction of data [1], some of which have been studied 
in detail, both according to theoretical models [2] and according 
to a series of experimental data [1], which are related, granted, 
to the technological parameters of the carrier rocket. 

The second method proposes solution of the problem of detection 
of the non-stationary states aboard the space vehicle. 

The solution on the selection of one of these two metnods 
may bo carried out with regard for the preliminary calculation of 
the coefficient of striction K, which characterizes the decrease in 
the number of data transmitted to the earth after carrying out of 
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the information sampling. 


For the first method, this magnitude is equal to 





) 


where N is the number of readings of the parameter with cyclic 
sampling during some time T, 


N is the number of readings remaining after the procedure of 

c 

elimination of the excessiveness, 


■ i>l is the coefficient which shows the increase in the volume 
of transmitted data, because of the necessity of transmission with 
striction of the supplementary function information (for recog- 
nition and/or temporary tie-in of the readings) . 


For the second method of sampling, the coefficient of striction 
is equal to 


K 


C ~ 



nc( 


s e c 
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where n is the number of non-stationary states of a given type de- 
tected during the time T, 


N is the length of the section, i.e., the number of readings 
s e e 

transmitted to earth for each detected non-stationary state. 


The number n includes some number of false responses, as a 
function of the detection algorithm. 

Examined in the present study is the first method of data 
sampling, and comparison and selection of data striction algo- 
rithms are carried out according to the results of experimental 
study of models cf striction devices. 

A separate study will be devoted to the task of searching for 
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a detection algorithm for realization of the second method. 

2. Procedure for Stu dying and Initial Data 

The study of the effectiveness of the data strictlon a „go _ 
rithms was carried out on the basis of the following three 
characteristics : 


1. The coefficient of strict ion K 

c 


/5 


2. The root-mean-square error in the readings after the regener- 
ation operation r - — n 


i • i 


3. The maximum error ir. the data after the regeneration oper- 
ation of the readings discarded as a result of the operation of 
the strict ion algorithm: Z. 


The first two character i st ics were determined as a result of 
calculations . Also selected were striction algorithms which give 
limited magnitude of the maximum error of regeneration. Keeping 
in mind the necessity of utilizing the results of the given analysis 
in the experiment on the study of shock waves in the interplanetary 
and near-planetary plasma, it was decided to carry out calculation 
of the characteristics of striction on some experimental data which 
describe the behavior of the plasma near the shock wave. Selected 
as such data in the first stage of the study were examples of 
measurements of the three components of the interplanetary magnetic 
field. These data were obtained using the SG-70 instrument on board 
the "Venera-9" space vehicle, during its flight in orbit as a 
satellite of Venus (November 1975-February 197b). Taken for 
analysis were sequences of communication with recording of the 
measurements in sections with a duration of from several minutes 
to several hours, with a frequency of sampling of each component 
of the magnetic field equal to 1 Hz. 
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According to tho procedure! qiven in [3], wo calculated the 
values? of the modulus- of the maqnet ie field, which were also 
utilized for subsequent analysis. 

In accordance with the posed problem, the selected data 
represented sections in which a near-planet shock wave is dis- 
tinctly observed, "calm" sections of the interplanetary medium 
up to tlu> front of the shock wave, and sections of "perturbed" 
plasma after the front of the shock wave. In addition, non- 
stutienar', areas were noted i r ' a number of sequences of communi- 
cation, in addition to the near-planet shock waves, which may be /i> 

interpreted as interplanetary shock waves or explosions in the 
plasma of differinq nature. 

The indicated frequency of samplinq makes it possible to not 
take into account the time constant of the maqnet ometer , which 
comprises about 0.2a sec. 

Shown in fiqures 1, 2, 3, and 4 (sequor A) are examples of 
a shock wave, "calm" and "perturbed" medium, and a potential 
interplanetary shock wave, and in fiqures 3 and u—— shock waves in 
sequences B and C. 

The data strict ion devices, which operate accordinq to 
diiterent alqorithms, were simulated on the computer in the form 
ot proqram modules. The data for proeessinq wore read out from 
maqnet ic tape, created as a result of the conduct of the staqes 
ot preliminary and primary proeessinq of the telemetric information, 
obtained in communicat ions sequences from the space vehicle. 

3 . Pat a S t r l c t i on A 1 qo r i t Jim s 

Selected for eornparison of effectiveness were three types of 
polynomial striction alqorithms, with an increasinq complexity of 
real izat ion . As is common knowlodqe, an increase in complexity 
'■loos not always lead to an increase in effectiveness, which is 
especially attributed to the striction of victual information; 
therefore, those studied were [1,2]: 


1. 7.ui*o order forecaster with floating aperture; 

2. zero order interpolator; 

3. first order interpolator. 

These algorithms ensure limitation of the maximum error of regener- 
ation by the prescribed values of *Z, where we will adopt Z as equal 
to 1 , or 2 r , keeping in mind that the interplanetary magnetic field 
usually is 3-10, ( , ~ 1 0 ~ oe) , the field beyond the shock wave is 

10-50,, and the accuracy of measurements of tne magnetometer is 0.5- 
1 , • 


We will examine the logical schemes of the striction algorithms 
and the methods of regeneration of the sequence of readings according 
to the condensed data. Insofar as there is no general ly-accepted /7 

terminology in the designation of algorithms, and because of the 
great number of their versions, the given logical schemes are 
sufficiently detailed. 


Zero Order Forecaster (PNP 

1. The reference reading X. is transmitted. 

2. Comparison of the reading difference modulus X. + ^ (initially, 

K 1, and the reference reading with a threshold Z is 

carried out. If |x. , -X. j <. Z actually, then we will move to 

i + k i 

step 3, or falsely — to step 4. 

3. The reading X , is considered surplus (passive), K=K+1 , we will 
return to step 2. 

4. The reading with the number i+k is considered significant (active), 
we will move to step 1, where i=i+k. 

Re gen e r at ion A lgor i thms 

The horizontal lines drawn to the right from each significant 
reading on the time scale of uncondensed sequence determine the values 
of all the intermediate surplus readings. The significant readings 
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which run in succession are connected by straight lines. 


Maxi-mini Zero Order In t erpolator (INP) 

1. The reading X. is taken. The value of Xj is assigned to the 
variables MAX, MIN, which designate the maximum and minimum values 

from the sequence of readings: 

MAX-X . 

1 

MIN=X. 

1 

2. The next reading X. + ^ is taken, where K-l initially. 

3. A search is conducted for a new maximum value in the group 

(MAX, X. .), and this value is assigned to the variable MAX. 

A search is conducted for a new minimum in the group (MIN, X [+[ ), 

and this value is assigned to the variable MIN. /8 

4. Verification of the activity in the group K+l of readings is 
carried out. If (MAX-MIN) 2Z actually, we will move to step 5, 
or falsely — to step 6. 

3. An evaluation is formulated of the group of passive readings, 
equal to 

KST = (MAX + MIN) / 2 . 

K is changed to K+l, and we return to step 2. 

6. The group of passive readings is terminated. The value of EST 
is transmitted, which represents the group of K-l passive 
readings, 1=1+K, and we return to step 1. 

Regeneration Algorithms 

Drawn on the time scale of uncondensed sequence of readings 
are horizontal straight lines, which determine the values of all 
passive readings to the left of each transmitted value prior to the 
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preceding transmitted value. 


First Order Interpolator (IPP) With T ransmis s ion of the Value 

1. The initial reference reading Xj is transmitted. We will move 
to step 3 . 

2. X i is transmitted as a reference reading. 


3. The upper and lower increments are calculated (with i=l, X.=X]): 

A - y, T + 7* - X; 

A i-hl + & > 

A 


B 


& 


n 


*-M “2 -X». 


4. The upper and lower values of the reading with the number i+k 
are formulated (initially, K=2): 


vB 

f h 


L-tK 




yH 

i+K = Xi + 


5. Verification of the activity of the reading X. + ^ is accomplished. 
If 


/9 


c T 


> k-., * & 

■ e fc-? fi 7 


actually, we will move to step 6, or falsely — to step 7. 


6. The reading X i + ] ^ is considered active, i=i+k-l. If K=2 , then 
X . =X , We will return to step 2. 

7. The reading X. . is considered passive. Verification of the 

. 1 k B 

mutual location of the upper value X. , and (X. ,+Z) is accorap- 

g 1 + N 1 K 

lished. If Xj + k'X. + j.+Z actually, we will move to step 8, or 
falsely — to step 9. 


8. Adopted as a new value of the upper value is 

* - ri v 

r 
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Tin 1 new value of the upper increment in calculated: 


d 





9. Verification of the mutual location of the lower value is 
accorupl is hod 

X Lr aml +k " ^ 

If actually, we will move to step 10, or falsely 

to step 11. 


10. Adopted as a new value of the lower value is 


A, 

vl 




(f . 


The new value of the ’ 'wor increment is calculated: 


11. The value of the last of the series K of surplus readinqs is 
cal cu luted 





K~K+1, we will return to step 4. 


Regener ati o n Algo ri thms /10 

The straight lines which connect the transmittable active 
readinqs on the time scale of uncondensed sequence determine the 
values of all the intermediate surplus readings. 


Also studied, in addition to these aluorithms, were two of 
their modifications, namely, the forecaster and interpolator of 
zero order, with a variable threshold. 


For PNP , the modification amounted to the following. With a 
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value of the reference reading X^IO,, no changes whatsoever took 
place. With X i 10 , , the magnitude 0.1 X i was utilized in step 2 
in place of the magnitude 2, i.e. / thu threshold was equal to 10V 
of the reference ading. 

Tor INP, the modification was somewhat more complex. In stop 
1, the assignment of the value of the initial reading X. to the 
value K8T is added. In step 4, if KSTslQ,, then the changes are 
not introduced, whereas if KST-10,, then the magnitude equal to 

q . i ( tUVX±M - IM ) utilized instead of the magnitude 7. This corres- 
> 

ponds to the designation of a threshold equal to 1 0 V. of the new 
value of the reading evaluation. 


Thus, with values of the parameter greater than 10,, the 
magnitude of the threshold 7 increases proportional to the increase 
in the values of the magnetic field, which leads to an increase in 
the coefficient of striction. It is natural that, in this case, 
the absolute value of the regeneration error will be greater the 
greater the magnitude of the reading. 


4. Results and Conclusions 

The values of the coefficients of striction, without taking 
into account the function information, i .e . , with ,~1 , and the 
root-mean-square errors of regeneration, obtained as a result of 
calculation according to three communications sequences from 
"Venera-9" for all of the studied algorithms, are given in tables 
1, 2, and 3. Recorded in each of these sequences is a shock wave, 
and the nature of the behavior of the magnetic field in "calm" 
sections is different. Sequences A and B were the most positive, /I 

and each of them contained long sections of "calm" and "perturbed" 
medium, with the "calm" section in sequence A having a much more 
variable nature than that in sequence B. Sequence B contained 
practically only small sections on both sides of the front of the 
shock wave, which had a substantial influence on the calculated 
values of the coefficients of striction. For the possibility of 
separate analysis, the coefficient of striction and the error were 
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calculated both according to the entire sequence, and accordin'* to 
the "calm" section, which includes the entire vo 1 nine of data on 
the interplanetary medium from the boundary of the sequence to the 
front of the near-planet shock wave, and the "perturbed" section, 
which contains data on the hot plasma from the front of the shock 
wave to the other boundary of the sequence. 

The calculations accordinq to the algorithms with a variable 
threshold were carried out only for sequence A, and only for the 
cane when, with a magnitude of the maqnetic field less than 10,, 
the threshold Z was equal to 1,. Comparison of the coefficients 
of striction shows that, of the three algorithms (FNP, INI’, IPP) , 
the qreatest coefficient of striction is provided by the zero 
order interpolator — -IdP. In sequence A, the maqnitude of the 

advantage, with respect to the first ordei interpolator , 
is equal to 1.43 with Z-l,, and 2 with 7 ,- 2 ., and, with respect 
to the zero order forecaster — 1.9 with Z 1 . , and 2.2 with 7 ,- 2 ,. 

The maqnitude of the advantaqe of the first order interpo- 
lator relative to the forecaster is small, which corroborates 
the conclusion, drawn earlier 12], on the inadvisability of the 
utilization of a sufficiently complex IPP algorithm. At the same 
time, the complication of the algorithm for the realization of 
the zero order interpolator, which is sliqht as compared with 
PNP, leads to a perceptible increase in the coefficient of striction, 
and makes it possible to show preference to it. 

Separate calculations of the coefficients of striction for 
"calm" and "perturbed" sections in sequences A an; 3 . D show that the 
effectiveness of striction is sharply different for them. In /12 

sequence A, in the "calm" section, the coefficient of striction 
increases by 1.9 times with Z-L,, and 6.7 times with 7 ,- 2 ,. For 
sequence B, the "calm" section is characterized by a practically 
constant average value and a small dispersion; therefore, the 
difference in the sections is more considerable for it: the 
coefficient of striction in the calm section increases by 2 times 
with Z = 1 ■ , and by 10 times with Z=2,. At the same time, the 
coefficients of striction throughout the sequence differ slightly. 
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The results for algorithms with a variable threshold show a 
slight (by 1.25 times for INP) increase in the coefficient of 
striction. This is explained by the fact that the basic contri- 
bution to the magnitude of the coefficient of striction throughout 
the sequence is made by the "calm" section, where the magnetic 
field is a; .eciably less than 10r# and the threshold Z does not 
change. This makes it possible to draw a conclusion on the inadvisa- 
bility of utilization of algorithms with a variable threshold in 
the given case. 

In addition to the magnitude of the coefficient of striction 
and the root-mean-square error, finer analysis of the distortions 
of the form of the signal, occurring after regeneration, is required 
for evaluation nj the effectiveness of the examined "quasiinvertible" 
methods for reducing the excessiveness. These distortions are 
illustrated by figures 7, 8, 9, and 10, where the same sections are 
shown as in figures 1, 2, 3, and 4, respectively, after regener- 
ation of the data, condensed by the zero order interpolator 
with a threshold Z=lj, and figures 11, 12, 13, and 14, for the 
threshold 2>. We will note some clear-cut features of the renener- 
ated data. First, the feature in the magnetic field, which looks 
like a potential interplanetary shock wave in figure 4, practically 
disappears in figures 10 and 14. In addition, that portion of the 
data directly preceding the sharp increase in the magnetic field 
in figure 1, i . e . , the front of the wave, is strongly distorted. 

In this connection, under actual conditions, it is inadvisable to /I 3 
make use, evidently, of the magnitude of the threshold Z = 2 Y - With 
a threshold Z=1 , , the magnitude of the coefficient of striction 
throughout the sequence did not exceed K c =6.38. T iking into account 
the fact that still another reading, which shows the current time 
of the measurement in the interval T, accompanies each reading of 
the parameter during realization of the algorithm of striction, we 
obtain i >2 in formula (1). In this case, the actual coefficient of 
striction for the sequence does not exceed the magnitude 3.19. The 
root-mean-square error of regeneration of the readings does not ex- 
ceed the magnitude 0.314, in this case, with a permissible maximum 
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error cf » 1 , . The obtained magnitude of the eoefficit : t of 
striction is not very groat, although one can nevertheless consider 
the presentation to the experimenter of information on the behavior 
of the parameter throughout the sequence, and not only in the area 
of the non-stationary state, an advantage of the use of polynomial 
algorithms . 

Evaluation of the competitiveness of polynomial algorithms 
of striction, relative to methods of detection, can be carried out, 
utilizing formula (2) , from the following considerations. With 
the transmission of data along a single loop of the planet satellite, 
the space vehicle twice intersects the near-planet shock wave, i . e . , 
n>2. Given a maximum volume of data, which may be recorded on board 
the space vehicle, and a volume of data which corresponds to a 
single section of the non-stationary area, having evoked response 
of the detection algorithm, we can determine the maximum number of 
recorded sections of the non-stationary area. Hence, we obtain 
the minimum value of the coefficient of data striction, which we 
can compare with the values which correspond to the polynomial 
algorithms. The comparison of the methods of striction will be 
carried out in the second part of the study. 

in conclusion, the authors express their gratitude to Candidate 
of Physical and Mathematical Sciences G. I. Zastenker for his 
useful discussions on posing of the problem and the results, to 
Candidate of Physical and Mathematical Sciences Ye. G. Yeroshenko 
for his submission of data and help in sampling of the actual / 14 _ 

information, and Candidate of Technical Sciences V. F. Babkin for 
his participation in the discussion of the results. 
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